Glass former units and transport in ion-conducting network glasses 
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A new theoretical approach is presented for relating structural information to transport properties 
in ion conducting network glasses. It relies on the consideration of the different types of glass forming 
units and the charges associated with them. Changes in the compositions of these units lead to a 
re-distribution of Coulomb traps for the mobile ions and to a subsequent change in long-range ionic 
mobilities. It is furthermore shown how measured changes of the unit compositions can be explained 
by thermodynamic modeling. The theories are tested against experiments on borophosphate glasses 
and yield good agreement with the measured data both for the compositional changes of the units 
and the variation of the activation energy. 
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The chemical composition of ion conducting glasses 
can be varied to a large extent and this offers many pos- 
sibilities to optimize these materials with respect to dif- 
ferent demands, in particular to high ionic conductivities 
1]. It is therefore important to get an understanding 
of the connection between the network forming structure 
and the long-range ionic transport properties. Consider- 
able progress has been made in the past to gain insight 
into near and medium range order properties of ion con- 
ducting glasses by various experimental probes such as 
X-ray and neutron scattering, infrared and Raman spec- 
troscopy, and solid-state NMR techniques (for a review, 
see (H). A challenge is to utilize this information for 
theoretical models of the ionic transport. One promis- 
ing route was suggested some time ago by building Re- 
verse Monte Carlo models of the glass structure based on 
diffraction data 3] and by further analyzingthese struc- 
tural models with the bond valence method [4| to explore 
the preferred diffusion pathways of the mobile ions. 

In this Letter we will present a new theoretical ap- 
proach, which is applicable to network forming glass 
structures and relies on the different network forming 
units (NFUs) that build up the host structure for the 
ionic motion (cf. Fig. [TJ. We argue that the charges as- 
sociated with the NFUs and the way how they are local- 
ized are of crucial relevance for characterizing the statis- 
tical properties of the energy landscape that govern the 
long-range ionic transport properties. To demonstrate 
the new approach we apply it to the mixed glass former 
effect in sodium borophosphate glasses, where detailed 
information on the NFU concentrations has been gained 
recently by MAS-NMR HQ, see Fig. H We first show 
how the observed changes of NFU concentrations with 
the borate-to-phosphate mixing ratio can be understood 
from a thermodynamic model. Then we will use this 
structural information on the NFUs to calculate changes 
of the conductivity activation energy upon the mixing 
ratio. 

In borophosphate glasses of composition j/Na20-(l — 



u)[xB203-(l — x)P20s] we distinguish seven NFUs as in 
[5j: the neutral trigonal B^ 3 ) units with three bridging 
oxygens (bOs) and zero non-bridging (nBOs), the nega- 
tively charged tetrahedral B^ 4 ) units with four bOs and 
zero nbOs, the trigonal B' 2 ' units with two bOs and one 
negatively charged nbO, and the tetrahedral phosphate 
units P", n = 0, . . . 3 with n bOs and (3 — n) nbOs, see 
the Fig. [TJ MAS-NMR measurements redrawn in Fig. [2] 
(symbols) show that, when starting the mixing from the 
phosphate rich side (x = 0), first the B^ 4 ) units replace 
P^ 2 ) units. This replacement continues until the B^ 4 ) 
concentration saturates at about x ~ 0.4. Above this 
mixing concentration the neutral B^ 3 ) units start to ap- 
pear, replacing now the neutral P^ 3 ^ units to keep the 
total amount of negative charge constant. This is needed 
to compensate the positive charge of the mobile sodium 
ions. With further increasing x, the behavior becomes 
more complex until at the boron rich side all NFUs are 
somehow involved in forming the network structure. 

To understand this behavior we developed a thermody- 
namic model, which is based on a hierarchy of formation 
enthalpies G(X) for the NFUs. The B< 4 ) are the most 
preferable NFUs for the charge compensation of the alkali 
ions, since they are most highly connected in the network 
with their four bOs. However, their concentration is lim- 
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FIG. 1: Sketch of the NFUs in borophosphate glasses. 
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FIG. 2: (color online). NFU concentrations in the glass 
0.4Na 2 O-0.6[a;B2O3-(l - x)P 2 O s ]. The symbols mark the 
MAS-NMR results from [jj and the open symbols at x = 1 
MAS-NMR measurements from The solid lines are the 
result of the thermodynamic model for A — > oo. 



ited, because the delocalized charge hinders them to come 
close to each other 0] ■ Defining in general the NFU con- 
centrations [X] as the fraction per network former cation 
(i.e. total number of units X divided by total number of 
B and P atoms) we set [B^] = min(x, [B( 4 )] sat ), where 
the saturation limit [B( 4 )] sat = 0.43 is chosen in agree- 
ment with the pioneering MAS-NMR results by Bray and 
coworkers Q and more recent findings [t| , as well as early 
theoretical modeling [l(| ■ For the other NFUs we choose 
the neutral P^ 3 ) unit as reference point, and introduce 
only one parameter A to describe the relative formation 
enthalpies 



SG(P {3) ) = SG(B^) = 8G(P {2) ) = 
<SG(P (1) ) = A, c5G(P (0) ) = 2A, 5G(BW 
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This choice expresses that the poorly connected and 
highly charged P^ 1 - 1 and P'°) units are increasingly less 
preferable compared to the better connected B^ 3 \ P^ 
and P*- 2 ) units and that the B^ 2 ) units are least likely 
due to their nbO and trigonal configuration, which makes 
it difficult to accommodate them within the network. In 
addition we need to take care of the total amount of neg- 
ative charge, being fixed by the sodium content due to 
charge neutrality, and the total amount of borate and 
phosphate given by x. This yields the constraints 



[B< 4 >] + [B< 2 >] 
[B^] + [B( 3 >] 

[p(3)] + [p(2)] 



[p(2) ] + 2 [p(i) 

[B< 2 >] = x 

[ P (D] + [P (0)] = (1 
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In a grand-canonical treatment we can assign the chem- 
ical potentials fj, q , fi B and /i P to these constraints (|2a|)- 
(Pc)l . respectively. Considering a set of sites to be occu- 
pied by the NFUs with mutual site exclusion, we obtain 
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|B< 2 >1 = 



exp[-fi x ] + 1 
1 



X = B (3) ;P (3) 
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(n)l 



exp[3A - /j, q - fi B ] + 1 
1 

exp[(2 - n)(A - fi q ) -fj, P ] + l 



(3a) 

(3b) 

n = 0,1,2 
(3c) 



where all energies are given in units of the thermal energy 
k^T and the chemical potentials have to be determined 
from Eqs. (|2lt [) -(|2c |) . 

Equations ([lj with the single parameter A describe 
the hierarchy between the formation enthalpies. Specific 
values for these enthalpies should be irrelevant as long as 
the system is in the low-temperature regime. To evaluate 
the behavior in this regime we solve the set of Eqs. (|2I3[) 
for A — > oo. The results shown as solid lines in Fig. m 
are in good agreement with the MAS-NMR data from 
ref. Q (diamonds), except for x = 1, where the mea- 
sured [B^ 4 )] is much smaller than the presumed satura- 
tion value [B( 4 )] sa t = 0.43, and correspondingly the [B( 2 )] 
value larger than the theoretical prediction. With respect 
to the deviation at x = 1, we note that MAS-NMR mea- 
surements reported by another group yield the data 
marked by the open symbols in Fig. [21 which are in bet- 
ter agreement with the theoretical predictions. On the 
basis of the thermodynamic model, one can, of course, 
reproduce the behavior found in ref. Q by assuming a 
lower saturation value [B( 4 )] sat for x — 1. Indeed, for the 
sodium borate glass, the maximal [B^ 4 ^] was found to be 
slightly smaller than 0.43 7], which can be explained by 
requiring that a bO cannot link two B^ 4 ^ units [lCf. How- 
ever, to describe all details, including different behaviors 
for different types of alkali ions, one needs to weaken 
this rule and allow for the formation of diborate groups 
12|. Let us note that by including such refinements it is 
also possible to model the [B^ 4 ^] max in the borophosphate 
system. To keep things simple we have focused on the 
essential idea and used the limit [B( 4 )] max ~ 0.43 here. 

Next we show how one can, based on the informa- 
tion on the NFU concentrations, successfully model long- 
range ionic transport properties. To this end we devel- 
oped a model, which we call the Network Unit Trapping 
(NUT) model. It relies on the following idea: the nbOs 
create localized Coulomb traps for the mobile ions, while 
delocalized charges, as those of the B^ 4 ) units, give a 
partial Coulomb contribution to several neighboring ion 
sites. In this way the structural energy landscape for the 
ionic pathways is modified with the mixing concentration 
x and this effect can be conjectured to govern the change 
of the activation energy E & (x) for the long-range ionic 
transport. 

To test this model we randomly distribute the NFUs 
with their concentrations from Eqs. ([3]) on the sites of 
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FIG. 3: Two-dimensional sketch of the NUT model. The 
arrows indicate charge transfer to ion sites as described in 
the text. 



a simple cubic lattice. These sites are called NFU sites. 
The mobile ions are considered to perform a hopping mo- 
tion between the centers of the lattice cells, which rep- 
resent the ion sites. An NFU a with k a > nbOs and 
charge (—z a e) adds a Coulomb contribution (—z a e/k a ) 
to k a randomly selected neighboring ion sites, as illus- 
trated in Fig. [3] Note that this implies that the dereal- 
ization of electrons belonging to the double bond in the 
charged P^ units is taken into account. For example, a 
p( 2 ) unit on an NFU site i induces a charge — e/2 at two 
randomly selected neighboring ion sites. The delocalized 
charge of a B^ 4 ) unit is spread equally among the neigh- 
boring ion sites, which amounts to set k = 8 for this unit. 
The neutral B^ 3 ) and P^ units give no Coulomb contri- 
bution. Finally, Gaussian fluctuations are added to the 
site energies in order to take into account the disorder in 
the glassy network [l3j | . In summary we can write for the 
energy of ion site i 



Ei = —En 



a, J 



where the sum over j runs over all neighboring NFU sites 
of ion site i. The occupation number is equal to one, 
if an NFU a on site j contributes a Coulomb contribution 
—z a e/k a to ion site i; otherwise it is equal to zero. The 
parameter Eq > sets the energy scale and the r)i are 
independent Gaussian random variables with zero mean 
and standard deviation a. Note that Eq is irrelevant as 
long we are interested in relative changes of the activation 
energy with x. Hence a is the only tunable parameter in 
the modeling. 

To determine the activation energy E a (x) we have cho- 
sen a lattice with 50 3 sites, occupied all NFU sites accord- 
ing to the occupation probabilities given by Eqs. ([3]), and 
the ion sites randomly with concentration y / ( 1 — y ) . Then 
Kinetic Monte-Carlo simulations with periodic boundary 
conditions and Metropolis transition rates [l4| were per- 
formed. After thermalization the time-dependent mean- 
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FIG. 4: (a) Arrhenius plot of the simulated Na + diffusion co- 
efficients D in 0.4Na 2 O-0.6[a;B2O3-(l - a;)P 2 5 ] for various 
x and a — 0.25. D is given in units of ua 2 , where v is the 
attempt frequency of the ion jumps and a is the lattice con- 
stant (mean jump distance). The slope of the regression lines 
yields the activation energies, (b) Comparison of the sim- 
ulated activation energy (open squares) with the measured 
conductivity activation energy from 0] (full squares). The 
solid line is a least square fit of a polynomial of sixth order to 
the open symbols. The open circle at x — 1 (connected with 
the dashed line) corresponds to the simulated E a , if the NFU 
concentrations from [f| are taken. 



square displacement R 2 (t) of the mobile ions and the dif- 
fusion coefficient D — lim^oo R 2 (t)/6t are determined. 
The diffusion coefficient is shown for a = 0.25 and vari- 
ous mixing concentrations in an Arrhenius plot in Fig. 2k,. 
From the slopes of the straight lines we calculated the 
activation energy E a (x), and the behavior of the nor- 
malized activation energy E a (x)]E a (0) is compared with 
the experimental results from [5J in Fig. The overall 
agreement between the theoretical (open symbols, solid 
line) and the experimental data (full symbols) is surpris- 
ingly good. Note that we needed to fit only one pa- 
rameter a to achieve this agreement. A significant dif- 
ference between the theoretical and experimental curve 
can be seen for x — > 1: while the theoretical E a (x) de- 
ceases monotonously with x, the experimental E a (x) fi- 
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nally rises for the sodium-borophosphate glass (x = 1). 
Interestingly, this rise is reproduced by the NUT model 
(dashed line), if instead of the NFU concentrations pre- 
dicted by Eqs. ([3]), the NFU concentrations measured in 
Q are used. In view of the discrepancies at x = 1 be- 
tween experiments discussed in connection with Fig. [21 
this calls for a reevaluation of the activation energy in 
the sodium borate system. 

In summary we have presented a new approach to re- 
late structural information to transport properties in ion- 
conducting network glasses. This approach is based on 
a consideration of the properties of the different NFUs 
building the network structure with respect to total 
charge and charge derealization. In addition we showed 
how MAS-NMR results for NFU concentrations can be 
understood from thermodynamic modeling. The poten- 
tial of our new approach is manifold, since one can apply 
it quite generally to other network glasses with different 
compositions. One immediate application, for example, 
could be the investigation of glass series with varying 
modifier content. It is known that the activation energy 
often shows a logarithmic decrease with the concentra- 
tion of mobile ions [l5| and it would be important to 
see whether this behavior can be captured by the NUT 
model. 
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